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Abstract 

We construct quantum effective action in spacetime with branes/boundaries. 

This construction is based on the reduction of the underlying Neumann type 
boundary value problem for the propagator of the theory to that of the much 
more manageable Dirichlet problem. In its turn, this reduction follows from 
the recently suggested Neumann-Dirichlet duality which we extend beyond the 
tree level approximation. In the one-loop approximation this duality suggests 
that the functional determinant of the differential operator subject to Neumann 
boundary conditions factorizes into the product of its Dirichlet counterpart and 
the functional determinant of a special operator on the brane - the inverse of 
the brane-to-brane propagator. As a byproduct of this relation we suggest a new 
method for surface terms of the heat kernel expansion. This method allows one 
to circumvent well-known difficulties in the heat kernel theory on manifolds with 
boundaries for a wide class of generalized Neumann boundary conditions. In par¬ 
ticular, we easily recover several lowest order surface terms in the case of Robin 
and oblique boundary onditions. We briefly discuss multi-loop applications of 
the suggested Dirichlet reduction and the prospects of constructing the universal 
background field method for systems with branes/boundaries, analogous to the 
Schwinger-DeWitt technique. 


1. Introduction 

Prospective goal of the present paper is to develop background field method [H |2] for 
brane models in gravity/string theory and cosmology. Current status of brane concept 
essentially relies on the analysis of quantum properties beyond the tree-level approx¬ 
imation. This is especially important in the problem of low strong coupling scale in 
brane induced gravity models which incorporate an infinite sequence of strongly cou¬ 
pled higher dimensional operators. Their calculation is most efficient in the language of 
quantum brane effective action. This, in its turn, requires application of the covariant 
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method of background field mq in which the propagators of the theory within certain 
approximation are calculable in the external (mean) helds of a generic form. 

The peculiarity of brane models is that their bulk propagators, rather than being 
dehned in inhnite spacetime with simple falloff conditions, are subject to nontrivial 
boundary conditions on branes. Since the helds are subject to dynamical quantum 
huctuations on timelike branes, these boundary conditions belong to the class of gen¬ 
eralized Neumann boundary conditions involving on branes the values of helds together 
with their normal and tangential derivatives. Finding such propagators (usually based 
on the method of images) is a very hard task, especially for helds with spins, when 
their Green’s functions have numerous spin-tensor indices. 

On the contrary. Green’s functions with Dirichlet boundary conditions are much 
easier to obtain - the method of images gives them as relatively simple linear com¬ 
binations of known propagators in spacetime without boundaries (dehned by mirror 
image continuation across the original boundary [3]). It turns out that Feynman dia¬ 
grammatic technique based on Neumann type Green’s function can be systematically 
reduced to that of the Dirichlet type, and the goal of the present paper is to develop 
the needed technique. This technique is based on the duality between the Dirichlet and 
Neumann boundary value problems, recently discovered in j3] at the tree level. Here 
it will be extended to the one-loop and multi-loop levels. 

The action of a (free held) brane model generally contains the bulk and the brane 
parts 

dX0(X)F(V)0(X) + J dx 

B b 

where the bulk ((i-l-l)-dimensional and the brane d-dimensional coordinates are labeled 
respectively by X = and x = x^, and the boundary values of bulk helds 0(X) on 
the brane/boundary are denoted by <p(x) 



- ip{x) K{d) ^{x) + i{x) cp(x) , ( 1 . 1 ) 


0(X) 


(f(x). 


( 1 . 2 ) 


The kernel of the bulk Lagrangian in given by the second order diherential operator 
F(V), whose derivatives X = dx are integrated by parts in such a way that they form 
bilinear combinations of hrst order derivatives acting on two diherent helds (this is 
denoted by F (V)). Integration by parts in the bulk gives nontrivial surface terms 
on the brane/boundary. In particular, this operation results for a symmetric operator 
F(V) in the Wronskian relation for generic test functions <pi^ 2 {X) 


j d''+^x(0iF(V)02-0iF(V)02) = - J PF02-01 IH^s) • (1-3) 

B dB 
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This relation can be regarded as a definition of the first order Wronskian operator 
W = hT(V) on the boundary/brane b = dB of spacetime bulk domain B. Arrows 
everywhere here indicate the direction of action of derivatives either on 0i or 02 - 

The brane part of the action contains as a kernel some local differential operator 
K{d), d = dx- Integration by parts here is irrelevant for our purposes, because b is 
considered to be either closed compact or having trivial vanishing boundary conditions 
at infinitely remote boundary db. j{x) plays the role of sources located on the boundary. 
The order of the operator K{d) in derivatives depends on the model in question. In 
the Randall-Sundrum model jB], for example, for certain gauges it is just an ultralocal 
multiplication operator generated by the tension term on the brane. In the Dvali- 
Gabadadze-Porrati model [ 7 ] this is a second order operator induced by the brane 
Einstein term, K{d) ~ □//! where /i is a very low DGP scale of the order of magnitude 
of the horizon scale, responsible for the cosmological acceleration jH]. In context of 
Born-Infeld action in D-brane string theory with vector gauge fields K{d) is a first 
order operator p. 

In all these cases the action (mi) with dynamical (not fixed) boundary conditions 
(p{x) for bulk fields naturally gives rise to generalized Neumann boundary conditions 
of the form 



(1.4) 


which involve normal derivatives of (p{X) contained in W (V) and generically also the 


tangential to the boundary derivatives contained in K{d) (and possibly in W (V)). 
These boundary conditions when imposed on all brane boundaries of the bulk along 
with regularity requirements at the bulk infinity uniquely define the Neumann type 
bulk propagator of the theory, brane-to-brane propagator, etc. and, therefore, uniquely 
specify all orders of perturbation theory for both tree-level and quantum effective 
action. 

Main result which we want to advocate here is that all the Neumann type ingredients 
of the perturbation theory can be systematically reduced to those of the Dirichlet 
boundary conditions 


(1.5) 


4>d — 0 . 


In particular, the tree-level brane effective action, obtained from (CH) by integrating 
out the bulk fields subject to boundary conditions reads as 



( 1 . 6 ) 


b b 


y) = - WGDW\\ix,y) + K{d) 6{x, y) 


brane 


(1.7) 
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with the brane-to-brane operator y) expressed in terms of the Dirichlet Green’s 

function GoiX, Y) of the operator -F(V) in the bulk. This expression implies that the 
kernel of the Dirichlet Green’s function is being acted upon both arguments by the 
Wronskian operators with a subsequent restriction to the brane. Double vertical bar 
indicates that both points of the operator kernel are restricted to the brane and labeled 
by corresponding low case letters. That is, if the embedding of the boundary/brane in 
the bulk is denoted hy X = e(x), then this explicitly means: 


WGDW\\{x,y) = W{Xx)Gd{X,Y)W{Xy) 


X=e{x),Y=e(y) 


( 1 , 8 ) 


It is obvious that F^^^^^{x,y) is essentially nonlocal, its local part being presented 
by the last term of (inD - the contribution from the brane. The Green’s function 
G brane ( 3:^5 u) of the brane operator. 


= 5{x,y), (1.9) 

is the brane-to-brane propagator of the bulk theory, and with the conventions of the 
above type this reads as the following expression for the brane restriction of the Neu¬ 
mann Green’s function Gx{X,Y) of F{S/) 



Gv|| i^X^y') GbraneI/)• (^’^0) 

The duality relations (inD and (fTTUl) were derived in j3] for a simplest case of K{d) = 
0. Below we generalize them to the case of a nontrivial K{d) and, moreover, extend 
them beyond the tree level. In particular, in the one-loop approximation we show 
that the functional determinant of the bulk operator T’(V) subject to the generalized 
Neumann boundary conditions dH factorizes into the product of the Dirichlet type 
determinant of T’(V) and the functional determinant of the brane-to-brane operator 
^brane boundary (d-dimensional) theory - the fact briefly reported in [^. This 

implies the following additive property for the one-loop effective action 

A-ioop = ^Tr^ Ini" = trluF^-'"-, (1.11) 

where Tr^ jv denotes functional traces of the bulk theory subject to Dirichlet and Neu¬ 
mann boundary conditions, while tr is a functional trace in the boundary d-dimensional 
theory. 

Gertainly, beyond tree level the effective action contains ultraviolet divergences, so 
that this one-loop Dirichlet-Neumann reduction property (HHD should be understood 
within certain regularization. Below we will use the dimensional regularization with 
the dimension d continued to the complex plane playing the role of regularization 
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parameter. In principle, other regularizations are possible, and their admissible types 
will be discussed in Sect. 2 below. 

As a byproduct of ()1.11|1 we suggest a new technique for surface terms of the local 
heat kernel expansion in spacetimes with boundaries. Heat kernel gives a proper time 
representation for the functional determinant of (pseudo)differential operators and, 
therefore, serves as a basic tool for the calculation of effective action in background 
held formalism D121 cni im ini In the presence of boundaries its local expansion is 
modihed by additional surface terms which amend easily calculable bulk terms well 
known in physics context as Schwinger-DeWitt coefficients. Calculation of these sur¬ 
face terms j2] presents a strong challenge of both technical and sometimes conceptual 
(nonperturbative) nature, especially for the so called oblique boundary conditions jT2] 
which contain derivatives tangential to the brane and arise, in particular, in Born-Infeld 
context j2]. Interestingly, Neumann-Dirichlet duality relations suggest an alternative 
method of their calculation, which in view of simplicity and universality has essential 
advantages as compared to the conventional approach of Uni d 1131121 • 

The paper is organized as follows. In Sect. 2 we derive the algorithms (inD, (innD- 
(dlD. In Sect. 3 we demonstrate them on several simple examples. In Sect. 4 we 
calculate several lowest order surface terms of the heat kernel for a wide class of gener¬ 
alized Neumann boundary conditions including, in particular, the case of oblique ones. 
In the concluding Sect.5 we briefly discuss the extension of the technique to multi-loop 
orders, its peculiarities in gauge theories and other problems of the curvature expansion 
in spacetimes with boundaries, which will be fully considered in the forthcoming papers 
uniini. In Appendix A the classical Feynman derivation of the gaussian functional 
integral in spacetime with boundaries is briefly revisited, and Appendix B is devoted 
to the derivation of the heat kernel for a special case of the brane-to-brane operator. 

2. Neumann vs Dirichlet problems 

We begin this section with specifying in more detail the structure of the second order 
bulk operator F’(V). As a kernel of the bulk action in (11.11) it should be symmetric 
and have the following general form 

F(V) = -dAa^^dB - b^dA + dA{b^f + c. (2.1) 

Its coefficients are some general coordinate dependent matrices acting in the vector 
space of (p{X) labeled by some spin-tensor indices which we do not specify here. With 
respect to these indices the coefficients and c are symmetric (a"^^)^ = = c. 

The Lagrangian of the bulk part of (HH) for this operator, containing the hrst-order 
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( 2 . 2 ) 


derivatives, is of the form 

(V)0= ^dA(pa^^dB(p- (pb^dA(p+ ^4>c(j). 

With one integration by parts, this Lagrangian differs by the total derivative term from 
the expression in which the operator i^(V) acts entirely to the right. For two different 
test functions this reads as 


</>! F 02 = 0i(F02) + -9a (01 w^02) 
in terms of the local Wronskian operator 
= a^^dB + b^ 

and can also be rewritten as a local Wronskian relation 

01 F(V)02 — 0iF(V) 02 = —(9a (01 W^4>2 — 01 hF^02^ • 
Integrating (Q over the bulk we have the equation 

J d'^+^X0iF02 = J d‘^+^X0i(F02) + J d'^(r0iW02 . 

B B b 


(2.3) 


(2.4) 


(2.5) 


( 2 . 6 ) 


It determines the boundary/brane Wronskian operator W which is given by the nor¬ 
mal projection of the local operator (E3D at the boundary (up to the measure factor 
involving the ratio of determinants of the bulk metric Gab and induced on the brane 
metric g^v), W= {^/g/\fG) W^- Similar integration of (12.hj) yields equation (II.dll of 
Introduction^. 

Now consider the functional integral in the brane model with the action dEH) 

^ =/9^0 exp(-F[0]), (2.7) 

where the integration runs over the bulk fields 0(W) and also over its boundary values 
ip{x), (III2D, on the timelike branes. Integration over the latter follows from the dy¬ 
namical nature of <p(a:) which are subject to independent quantum fluctuations. This 
gaussian path integral equals 

Z = (DetGiv)'''"exp(-F[0iv]), (2.8) 

^Note that the Wronskian relation specifies W(V) only up to arbitrary symmetric operator 
acting the boundary (like n{d)), while Eg. (12.till fixes it uniquely. 
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where 07 v is a stationary point of the action (HU satisfying the following problem with 
the inhomogeneous Neumann boundary conditions 


F(V) MX) = 0, 

(ly +k)(/)7V +j{x) = 0, 


(2.9) 


and Gat is the Neumann Green’s function of the bulk operator - the solution of the 
following problem 


F(V)G^(X,F) = 5(X,F), 
{W +K)Gn{X,Y) =0. 


( 2 . 10 ) 


For completeness in Appendix A we present the derivation of the gaussian integral 
(EHl) by Feynman’s method UBI which clearly shows how the boundary conditions 
enter the calculation of the preexponential part of this algorithm. This derivation, in 
particular, gives the variational dehnition of the corresponding functional determinant 
which goes far beyond its matrix (hnite-dimensional) analogue. It is important that the 
boundary value problem (El naturally follows from the action (El and Wronskian 
relations for F’(V), because the variation of the action is given by the sum of bulk and 
brane terms 


5S[(p] = / dX dcj) {F(j)) + dx 5^p {W(p + Kcj) + j) 


' B 


( 2 , 11 ) 


which separately should vanish (remember that the action should be stationary also 
with respect to arbitrary variations of the boundary helds 5^). 

The solution of (El has the following form in terms of the Neumann Green’s 
function 


0jv(X) = - dyGN{X,y)j{y) = -GnU, 


( 2 . 12 ) 


GN{X,y) = GN{X,Y) 


Y=e{y) 


and the stationary action as a functional of the boundary source j{x) equals 


-S') ] = ^ / dX 0 (F 0) + I dx( +k) (p + jcp), 


IB 


= -- dxdyj{x)GN{x,y)j{y) = -- JGnW j 


X=e{x),Y=e(y) 


= G 


N\ 


(2.13) 

(2.14) 


GK(x,y) = G,,(X,Y) 


7 







Here to simplify the formalism we used condensed notations by omitting the sign of 
integration over boundary/brane coordinates'^. Thus finally we have 

Z = (DetG7v)^''%xp JGnW j)- (2.15) 

Alternatively one can calculate the same integral by splitting the integration pro¬ 
cedure into two steps - first integrating over bulk helds with hxed boundary values 
followed by the integration over the latter 

j D4> {...) = j I Dct> (2.16) 

4>\ = v 

This allows one to rewrite the same result in the form 
^ = / d(pZ{(p), 

Z((fi)= J D(p exp{-S[(j)]), 

= 

where the inner integral over bulk helds in view of gaussianity is again given by the 
contribution of a saddle point (j)o 

Z{ip) = (DetGD)^'^^exp(-S'[0D]). (2.19) 

This saddle point conhguration satishes the problem with the inhomogeneous Dirichlet 
boundary conditions 

F(V) 0d(X) = O, 

0z)|=V9(a:), (2.20) 

and the preexponential factor of is given by the functional determinant of the 

Dirichlet Green’s function subject to 


(2.17) 

(2.18) 


F(V)Gd{.X,Y) = S(X,Y), 
Gd(X,Y) =0. 


( 2 . 21 ) 


^We will never use this rule for bulk integration which will always be explicitly indicated together 
with the corresponding integration measure. It is useful to apply this rule for integral operations 
on the brane, though, because these operations never lead to surface terms and in our context have 
properties of formal matrix contraction and multiplication. In the case of the one-dimensional bulk 
this rule applies literally, and it can be extended to higher dimensions without any risk of confusion. 




In terms of this Green’s function and using condensed notations we have 


MX) 


dy Gd{X,Y)W 


ipiy) = -GdW\ 

Y=e{y) 


( 2 . 22 ) 


^[4>d] = -^ dxdy(p{x) - WGdW {x,y) + Hi{x,y) (p{y) + Ji^) M) 



- WGdWW + K 




(2.23) 


where W GdW \ \ is dehned by Eq. m in Introduction. Note that the last expression 
is exactly the tree-level brane effective action obtained from the original action 
by integrating out the bulk helds subject to boundary conditions (p{x), 


^ brane [ ^ ] ^ [ 4^D [^] ] • 


(2.24) 


Substituting with into (ITT7I) we again obtain the gaussian integral 

which is saturated by the saddle point ipo of the above brane action (E?l) 




- WGdW\\ + k 


1 -1 


(2.25) 


and the hnal result reads 

Z=fDetGz)j ( det[-W GdW \\ +k]] 

X exp j [-IEGdIEII + , (2.26) 

where det denotes the functional determinants in the d-dimensional boundary theory. 

Comparison of its tree-level and one-loop (preexponential) parts with those of (j2.15|) 
then immediately yields two relations 

Gn\\ = [—WGdW\\+k] ^ = G brane) (2.27) 

fDetGivj = f DetGDj (det[- WGdW\\ + k]) . (2.28) 

They underly the algorithms (USD - (EZI), (fTTUD - lfrTTD advocated in Introduction and 
give the possibility in a systematic way to express all the quantities in brane theory in 
terms of the objects subject to Dirichlet boundary conditions. 

As mentioned in Introduction, the relation between functional determinants 
should be understood within some ultraviolet regularization. It should regulate the 
both bulk functional integrals (j2.7|) - (j2.8|l and (j2.18|) - (j2.19|l as well as the boundary 
functional integral over ip = (p{x) (of d-dimensional theory) in (I2.16j) - (I2.17D . The 
simplest procedure which regularizes all three integrations without violating the basic 
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relation (innD consists in the analytic continuation in d to the domain of convergence 
of relevant Feynman integrals. Other types of regularization do not seem to violate 
(irrai too, however they can qualitatively change the setting of the boundary value 
problem underlying the obtained algorithms. For example, the regularization by higher 
derivatives, as well as certain versions of the Pauli-Villars regularizations, increases the 
order of differential operators. This changes the order of the normal derivative in the 
boundary conditions (El and even the number of the latter, so that the Dirichlet- 
Neumann reduction stops working or has to be essentially modihed. For this reason, in 
what follows we will use dimensional regularization as the simplest and most efficient 
scheme. As we shall see in Sect. 4, it correctly generates the ultraviolet hnite heat 
kernel underlying the calculation of functional determinants by the Schwinger proper 
time method and, thus, conhrms the validity of the chosen regularization technique. 


3. Simple examples 


3.1. One-dimensional problem 

The simplest case of the Dirichlet-Neumann duality relation can be demonstrated on 
the example of the one-dimensional Sturm-Liouville problem on the segment of hnite 
length = I 


F(V) = 




y-<y< y+- 


Its Wronskian operator on the two boundaries of this segment is given by 


W 


= 


and the Dirichlet and Neumann Green’s functions are correspondingly 
sinkm(y'-y+) sinkm(y-y_) 

Gniy, y ) =- r-r — j - 0{y -y) + [y ^ y), 

m sinhm/ 

Gwl!/, y ) =- z-—rz7, - S(9 -a) + (a «j/). 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


m sinhm/ 

Since the boundary of the one-dimensional bulk consists of two points y±, the full 
brane operator Ha has the form of the two-dimensional matrix with the elements 
corresponding to the ± entries on the two zero-dimensional ’’branes” (for simplicity we 
take the case of k = 0 ) 


T brane 


m 


sinh ml 


coshm/ —1 
—1 coshm/ 


(3,5) 
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On the other hand, the restriction of the Neumann Green’s function (Q to the bound¬ 
ary is given by 


Gn 11 — 


m sinhm/ 


cosh ml 1 
1 cosh ml 


(3.6) 


This is a matter of a simple verihcation to check that these two matrices are inverse 
to one another which is just the relation (HH). 

To check the one-loop duality relation one can write the Dirichlet and Neumann 
functional determinants of F(V) as products of eigenvalues of the corresponding spec¬ 
tra. Interestingly, for this simple problem the Dirichlet spectrum 


F(V).(.f = 4>°(y±)^o, 

(3,7) 

7r2 P 

= ^ + fc = l,2,3,... 

(3.8) 

coincides with the Neumann spectrum 

= dyct>Uy±) = ^. 

(3.9) 

Af - ,2 /c - 0,1,2,3,... 

(3.10) 


except one eigenmode k = 0. This constant mode, {y) = const, is just absent in the 
spectrum of the Dirichlet problem. Therefore 


Detjv T’(V) —Det£)T(V)Ao, (3.11) 

and this immediately conhrms the relation in view of the fact that 

detF^''"“ = m2 = Ao. (3.12) 


3.2. Half-space with Killing symmetry in extra dimension 

Another example is of held-theoretic nature. It corresponds to {d + l)-dimensional 
half-space with the d-dimensional boundary plane. Let the operator be given by the 
(d -|- l)-dimensional d’Alembertian with mass and let the boundary be located at the 
position 1 / = 0 of the extra-dimensional coordinate y = 

F(V) = =m‘^ -D-d^ □ = 

= (x^, y), y>0, 

X\b = (a;,0) 
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(3.13) 

(3.14) 

(3.15) 






The h-dimensional part of the full d’Alembertian can in principle be curved and non- 
trivially depending on x-coordinates. We only assume the possibility of separation of 
variables, so that y is a. Killing direction in the bulk. 

For such a setting the Wronskian operator is given by the normal derivative with 


respect to the outward-pointing normal and equals 

W= -dy, (3.16) 

while the exact Dirichlet and Neumann Green’s functions are 

Gd.k(!/. 9') = i T . (3.17) 

A = (3.18) 

where the minus and plus signs refer respectively to the Dirichlet and Neumann cases. 
According to (HID the brane-to-brane operator equals 

^brane ^ (?/, 2/') ^ U=.'=0 = A, (3.19) 

while from (ITT7I) 

GwII =G)v(0,0) = ^, (3.20) 

which conhrms the relations (HID and (ITTTUl) . 


To check (HUD let us write a variational dehnition of the functional determinants 
for both Dirichlet and Neumann boundary conditions with respect to general variations 
of the d-dimensional part of the full operator 6F = —6 □. We have 

poo 

6 \n Dei Gd,n =Gd,n SF = ti / dy GD,Niy,y) SO, (3-21) 

Jo 

where we decomposed the {d + l)-dimensional functional trace into the operation of 
integrating over y the coincidence limit of the corresponding ^/-dependent kernel and 
the d-dimensional functional trace tr. Then, substituting (IXT7I) we have 

/•OO T 

(5 (In Det Gjv — In Det Go) = tr / d?/— <5 □ = —<5 In det A, (3.22) 

Jo A 

which, in view of dsm, in the inhnitesimal variational form fully conhrms (HUD- 
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4. Boundary terms of the local heat kernel expan¬ 
sion 

Neumann-Dirichlet duality can be used for the calculation of the boundary terms in 
the local expansion of the heat kernel. In the absence of boundaries with trivial falloff 
conditions at inhnity, the heat kernel 

K{s\x,y) = ( 4 . 1 ) 

turns out to be a very efficient tool of the covariant diagrammatic technique for quan¬ 
tum effective action in curved spacetime and in external helds of a very generic form 
HI cni HB US- Its efficiency is based on the possibility of expanding this kernel in 
asymptotic series in integer powers of s —*■ 0 with the coefficients a„(X, Y) which sat¬ 
isfy simple recurrent equations. These coefficients, often called in the physics context 
the Schwinger-DeWitt or HAMIDEW coefficients, can be explicitly found in the co¬ 
incidence limit y = X as local invariants built in terms of spacetime curvature of the 
bulk metric Gab, fibre bundle connection and other background helds. Thus, they give 
rise to local low energy expansion of the effective action in inverse powers of the mass 
parameter > oo, when the inverse propagator of the theory is supplied with the 
mass term, E(V) — >■ E(V) -|- 

In spacetime with boundaries the situation becomes more complicated, because, 
similarly to Green’s functions, the heat kernel should be obtained from that of the 
inhnite spacetime by the method of images. This leads to the expansion of the func¬ 
tional trace of the heat kernel in half-integer powers of the proper time parameter 

jSl ca 1201112 

CXD 

(.">!„ + ) , ( 4 . 2 ) 

Together with the volume (bulk) terms of the Schwinger-DeWitt type^, 

An = J d‘^^^XG^/^{X)tTan{X,X), ( 4 . 3 ) 

M 

this expansion acquires surface integrals at the boundary Bn /2 which are built of local 
invariants incorporating also such local characteristics of the surface as its extrinsic 
curvature and the curvature of the induced metric 

Bn/2= [ d‘^xg^^‘^{x)bn/2ix). ( 4 . 4 ) 

_ dM _ 

^In this section we denote the functional trace in the (c? -I- l)-dimensional bulk and on the d- 
dimensional boundary respectively by Troj^+i) and Tr(£;), while the notation tr is reserved for the 
trace over spin-tensor indices of matrices. The latter are denoted by hats. 
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Here G{X) and g{x) denote the determinants of the bulk and brane metrics, so that 
tra„(X, X) and bn/ 2 {x) turn out to be the bulk and boundary scalars. 

In contrast to the bulk Schwinger-DeWitt coefficients a„(X, X) which are univer¬ 
sal and independent of the type of boundary conditions, the surface terms essentially 
depend on the latter, and their calculation is much less universal and often very cum¬ 
bersome. For the operator of the form 

F( V) = - VaVb -P+IrI (4.5) 

D 

with second order covariant derivatives forming a covariant {d + 1)-dimensional 
D’Alembertian in the metric Gab and P denoting some matrix-valued potential term, 
few lowest Schwinger-DeWitt coefficients read as 

ao(X,X) = i, ai(X,X)=P, .... (4.6) 


The corresponding surface terms for the Dirichlet and simple Neumann boundary con¬ 
ditions 

b^'^{x)=0, (4.7) 

bi/ 2 {x) = tr i, ^ tr i, (4.8) 

&f’^(x) = ^Xtri, (4.9) 


involve the trace of the extrinsic curvature of the boundary K = 01121 . 

For the generalized Neumann (Robin) boundary conditions 


(V„-F)0(X) 


dM 


= 0 


(4.10) 


the last of the above coefficients is modihed by the matrix-valued potential term S 


b^{x) = tr 


2S+-KI 

3 


(4.11) 


This modihcation is even more sophisticated in the case of the so-called oblique 
boundary conditions OBI, which include the tangential to the boundary (d-dimensional) 
covariant derivatives D,, 


Wn-P^D,--{D,P^)-S)ct>{X) 


= 0 . 


dM 


(4.12) 


These derivatives enter the boundary conditions with the dimensionless matrix-valued 
vector coefficients P^. For the generic case their contribution to bn /2 is not known, but 
in the case of commuting matrices 


'P'^,P‘'] = 0 


(4.13) 
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lengthy calculations of [r]i[ni[El lead to the following expressions for few lowest-order 
surface densities (see also nmg for higher order feri/ 2 ) 


TT 


bi/2ix) = ^tr 
fef (x) = tr 


VTTr2 


- 1 


-^^S+-K1 
L 1 + r2 3 


1 + r2 


arctanh\/— 


r r 

p2 


where 


n = r^r,, 


(4.14) 

(4.16) 

( 4 . 16 ) 


It is important to note that these matrix functions are nonpolynomial in because 
of the dimensionless nature of these matrices. This means that their contribution to 
any given surface coefficient bn /2 cannot be obtained by the perturbation theory in F^, 
which is the main reason of difficulties in their derivation. 

Let us now give a simple derivation of these coefficients by the technique of the 
previous section. To begin with, note that the Wronskian operator for (031) is defined 
by the normal derivative with respect to outward-pointing normal to the boundary 


W= V„. V„ = 


(4.17) 


Introduce the covariant operator k{D) corresponding to the generalized Neumann 
boundary condition (imi) and the d-dimensional brane action, from which these bound¬ 
ary condition can be obtained by the variational procedure of Sect.2. Obviously, they 
read 

K = k{D) = -r>^D^ - 1 {D^fn - S, (4.18) 

^ + S)(p{x), (4.19) 

provided the following symmetry property of matrices F^ and S hold (which we assume 
in what follows) 

^ ^ ^4 20 ) 

According to the Shwinger-DeWitt proper-time method cmni the functional deter¬ 
minants of massive operators are given by the proper time integrals of the corresponding 
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heat kernels. In view of the representation (|12D this gives the following inverse mass 
expansions in the Dirichlet and generalized Neumann cases 


Tr^i^Arln 


F(V) + 



g—sF(V)—s m? 


/ m? Y'^+i)/2 
V dvr / 


n=0 


r (n 




m 


2n 


An+^ 


r(n^) 


B 


n=l 


m'' 


D,N 

nl2 


( 4 . 21 ) 


The heat kernel trace ( 1121 ) is always ultraviolet hnite, and the one-loop divergences 
originate from the proper-time integration diverging at s = 0. In dimensional regular¬ 
ization (with d analytically continued to the complex plane) they arise here as gamma 
function poles for first few terms of nonnegative powers in the mass parameter (in even 
dimension d -|- 1 for bulk Schwinger-DeWitt coefficients and both even and odd d -|- 1 
for surface terms). 

The Dirichlet and Neumann expressions (ICTl) differ only by the contributions of 
surface integrals, their total difference, on the other hand, being defined from the 
duality relation ma. Subtracting the Dirichlet version of (j4.2Hl from the Neumann 
one, one therefore obtains 


Tr(rf) In F 


brane 


1 .^y/2 “ 

2 V dvr / 

n=l 


hN 


hD 

%/2 


TT 


(d.22) 


where the brane-to-brane operator is defined by Eq. dnj in Introduction for a 

particular case of the boundary operator (iroi) . Thus, the difference of the boundary 
terms for Neumann and Dirichlet cases can be disentangled from the functional deter¬ 
minant of iT't’i'ane relevant orders of the 1/m-expansion. In view of the ultraviolet 
hniteness of the heat kernel (finiteness of the coefficients the divergences of this 

functional determinant should have the structure of gamma-function coefficients in the 
right hand side of which serves as a consistency check of the whole procedure. 

Though the operator is itself not exactly known, all we need is its inverse- 

mass expansion which is equivalent to the bulk-brane curvature expansion of the Dirich¬ 
let Green’s function of F(V). Moreover, it is defined on the brane manifold without 
a boundary (or with trivial regularity conditions at inhnity). All this essentially facili¬ 
tates the solution of the problem. As we show in the following examples, choosing the 
brane operator ()1.7j) for k{D) of the form (jd.lhj) we easily calculate the lowest order 
surface terms, and this procedure can undoubtedly be extended to all bn/ 2 - 
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4.1. Simple Neumann boundary conditions 

For simplest Neumann boundary conditions with k{D) = 0 the brane-to-brane operator 
(HU was obtained in the leading order approximation in Sect.3.2, Eq. dnm. It reads 

^brane = ^ ^ ^^2 - Q + 0[R, K], (4.23) 

where 0[R,K] denotes corrections due to the bulk curvature and extrinsic curvature of 
the boundary. Therefore, the inverse mass expansion for its determinant is dominated 
by 

Tr(,) InF^’’-'^- = ^ Tr(,) -D)+0[R,K] 

1 / \ dl‘1 f ^ 

= -2\^) n-dm j d^xg^/\il + 0[R,K] (4.24) 

with 0[R,K] = From the n = 1 term of (j4.22jl it follows then that 

^1/2 - ^f/2 = b (4.25) 

which fully agrees with (USD. The dependence of (j4.24j) on dimensionality (yielding the 
logarithmic divergence for even d) is exactly the same as in the n = 1 term of (j4.22jl 
which, as expected, guarantees the ultraviolet hniteness of the obtained difference 
^ 1/2 ^ 1 / 2 - 

4.2. Robin boundary conditions 

For Robin boundary conditions ()4.10|) with k,{D) = —S 

= Vm‘^ - n - S + 0[R, K], (4.26) 

The functional determinant of this operator can be obtained by perturbation theory in 
the dimensional quantity S 

IV,Inf = i Tr,j, ln(m" - □) - Tr,,, + 0\R,K,S^]. (4.27) 

Only the hrst order in S contributes to in the Robin case, and to zeroth order in 
the curvature this term equals 
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Therefore, from the n = 2 term of it follows that 

-h^ = 2iiS + 0[K], (4.29) 

which fully agrees with (031) and (mm). 


4.3. Oblique boundary conditions 

For simplicity consider the case of oblique boundary conditions (iroi) with S' = 0. 
Then the surface operator (iraii is k{D) = + 0[Dr ], and the brane-to-brane 

operator reads 

^brane _ ^^2 _ Q _ + O [ i?, JF, DT ], (4.30) 


where together with curvature terms we disregard terms with covariant derivatives of 
With the same accuracy 







0 


y=x 


(4.31) 


Here we denote the proper time parameter by r to emphasize that it has the dimension¬ 
ality different from s (length rather than length squared), which in its turn is explained 
by the dimensionality of (ICTl) . 

In contrast to the Robin case further calculations cannot be performed by pertur¬ 
bations in powers of the term because the vector coefficient F'^ is dimensionless, 

and the perturbation theory in this term will not generate asymptotic expansion in 
inverse mass^. Instead, with the same accuracy of zeroth order in the curvature one 
can disentangle this term in the exponential as a matrix-valued shift operator 


exp 


-rVm^ - □ + V, 


5{x,y) 


y=x 


= e 




K{T\x,y) =K(T\x + Tr,x) 


(4.32) 


y=x 


acting on the heat kernel of the operator \/m? — □, 


K{T\x,y) = exp 


—r — □ 


S{x,y). 


(4.33) 


^In other words, the leading symbol of the pseudodifferential operator - the highest (first) 

order in derivatives part of jt'_ jg given by \J and it should be treated as a whole 
without breaking into pieces. 
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Note that in view of the commutativity assumption Km this nontrivial matrix-valued 
function is uniquely dehned without any matrix ordering prescription®. 

As shown in Appendix B, in flat {d -|- l)-dimensional space without boundaries 

K{T\x,y) = 2T(^^j K(^d+i)/ 2 {Z), (4.34) 

where 


Z = Z{T\x,y) = m\/\ X — yl"^ + 


(4.36) 


and Ku{Z) is the modihed Bessel (MacDonald) function. Because of 


Z{t\x + Tr,x) = mr y 1 + r'^ 
the r-integration in ()4.31j) gives 


(4.36) 


Trill 



OO 



0 


K{t\x + tF, x) 


/ m? \(‘^+^h2 2 
V 27r / m 




K{d+i)/2{Z) 

Z{d+l)l2 ’ 


whence in view of the integral ( Eq. 6.561.16 of pT] i 


(4.37) 


OO 

j dx x-'^K.ix) = 2-"-'v^ r (4.38) 

0 

we hnally have 

Tr,„lnF'>'“' = -i(^)''''T(-l) j (4.39) 

Therefore, from the n = 1 term of Km it follows that 

b ?/2 - bf /2 = 0r tr ^ , (4.40) 

Vl + F^ 

which fully agrees with (imi) for oblique boundary conditions with ^ = 0. Similarly, 
one can check the S'-dependent term of for the case of nonvanishing S. 

The approximation of zero bulk and brane curvatures in all the above examples can 
be used as a starting point of the perturbation theory in 0[R,K] (and in other di¬ 
mensional background held quantities like D^F^ and P). Then the Neumann-Dirichlet 
duality method of the above type will give extrinsic curvature terms of (h:t^ and all 
higher order surface terms in the heat kernel expansion (Q. 

^Apparently, the commutativity assumption can be removed. For non-commuting matrices equa¬ 
tion will still be valid under the symmetric matrix ordering prescription in the right hand side. 

This symmetrization follows from the symmetry of the Taylor series coefficients. 
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5. Conclusions 


( 6 , 1 ) 


Thus, the algorithms (HHl), (fTTUli - ifrTTl) allow one to reduce calculations of brane 
effective action to those of the Dirichlet boundary conditions. This reduction technique 
can obviously be extended to multi-loop orders by applying the same trick of splitting 
the functional integration into two steps, as in - (jnHD, also in the nonlinear 

case. The resulting Feynman diagrammatic technique from combinatorical viewpoint 
is not so simple as in (ll.ll|l . but still manageable. In addition to the bulk Dirichlet 
type propagator Gd it has the brane-to-brane propagator G brane ^ the Green’s function 
of iT't’rane modihcation of Feynman diagrams, therefore, consists in the insertions 
into the bulk diagrams of the Dirichlet type the lines connecting bulk vertices to the 
brane by the brane-to-boundary propagators 

SMX) r (xY\u, 

r , X = -Gd(a, y) w 

d^[y) Y=e{y) 

(cf. Ea. fl2.22|l l and also developing the d-dimensional diagrammatic technique on the 
brane with the propagator Gbrane- This technique will be considered in more detail 
in cni- In gravitational brane models the question of gauge invariance (especially 
with respect to general coordinate transformations) becomes very important, while 
at present even the details of Faddeev-Popov gauge hxing procedure in spacetimes 
with boundaries are not clearly studied [5]. Therefore, gauge properties of Neumann- 
Dirichlet duality will be a major aspect of ^H], where gravitational Ward identities in 
brane models will be established (they are briefly reported in [^). 

As a byproduct, the Neumann-Dirichlet reduction technique suggests also a new 
method of calculating surface terms of the heat kernel expansion. Given the Dirichlet 
type terms, those of the generalized Neumann case can be obtained from (j021)- This 
allows one to circumvent the limitations of the conventional method for these terms. In 
particular, we were able to recover correct expressions for few lowest surface contribu¬ 
tions to the Schwinger-DeWitt expansion and, in case of oblique boundary conditions, 
perhaps even extend their validity beyond the commutative case (see footnote 

after Eq. on the symmetric matrix-ordering prescription). 

Of course, the success of Neumann-Dirichlet reduction programme depends on our 
ability to hnd the Dirichlet Green’s function Gd and the corresponding brane-to-brane 
operator xhe latter is a nonlocal pseudodifferential operator, so the problem 

of efficiently handling its nonlocality arises. For a wide class of problems F 
found in the zero curvature approximation as 


brane 


was 


1 brane 


(D) = Vm^-n + K(D) + 0[R,K]. 


(5.2) 


Despite nonlocality, local expansion of its functional determinant is still manageable. 
For ultralocal k{D) = —S, as in Robin case, it is easily available by perturbations. It is 
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more complicated in the case of oblique boundary conditions, when the leading symbol 
of gives rise to the propagation off light cone on the brane (in the physical 

theory with the Lorentzian signature) - the phenomenon called generalized causality 
in cni'- Finally, in brane induced gravity models with the operator k{D) r\j 
generated by the brane Einstein term (/i is the DGP scale [3 Ej), the calculational 
technique still has to be worked out QH. 

Final comment of this paper concerns the (bulk and boundary) curvature expan¬ 
sion of the Dirichlet Green’s function. Glassical method of images for this expansion 
is known jS], and this method for the Dirichlet case is much simpler than for the gen¬ 
eralized Neumann boundary conditions. Still further efforts are necessary to convert 
it into a regular systematic calculational scheme comparable in its universality to the 
Schwinger-DeWitt technique in spacetime without branes/boundaries |l()l ITT] . The 
progress in this direction will be reported in forthcoming papers naim, and here it 
remains to express a hope that at least the general shape of background-field method 
for quantum effective action in brane theory becomes visible. 

A. Gaussian path integral in spacetime with bound¬ 
aries 

Feynman’s calculation JH] of the gaussian functional integral (EHl) is based on the 
integral 

Z[F, J] = J D<j) exp (-F[0, J]). (A.l) 

flere instead of the source j (x) at the boundary, as in the action has the source 

J{X) to the integration field 0(x) in the bulk 

S14>,J]^ jjx (t 4,{X)F(X)4>{X) + J{X) 4,{X) 

-I-- / dx(p{x) K{d) ip{x). (A.2) 

2 JdB 

To hnd the dependence of (E3 on J{X) consider the stationary point of this action 
with respect to arbitrary variations of (p{X) both in the bulk and on the boundary. 

^Related phenomenon of loss of strong ellipticity, which is the unboundedness of the operator 
spectrum from below in the Euclidean theory [22] or the presence of ghost modes in the Lorentzian 
case, takes place when 1 + acquires zero or negative eigenvalues leading to singularities in the 
algorithms (ITTni - (imil 
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Similarly to (iznD this field satisfies the generalized Neumann boundary value problem 
F(V) MX) + J{X) = 0, 


{W +k)(Pn 


= 0 . 


Now make the shift of the integration variable in (El) by 07V 
0 = 0w + A. 


(A.3) 


(A.4) 


Under this replacement the action decomposes in the part S'[ A, 0] quadratic in A and 
the part independent of A. Linear in A term is absent (both in the bulk and on the 
boundary) in view of the stationarity of the action at 0Ar, so that 


^[0,J] = A[A,O] + ^[07v,J], 

A[0jv,J] = ^ [ dXdY J{X)Gn{X,Y)J{Y). 
2 Jb 

Therefore 

Z[F,J] = Z[F,O]exp(-A[07v,J]), 


(A.5) 

(A.6) 

(A.7) 


which justifies the exponential (tree-level) part of ()2.8|1 . 

To find the prefactor, consider the variation of the integral (El) at J = 0 with 
respect to the operator -U(V) and make the following set of obvious identical transfor¬ 
mations using the above equations (lA.bfl and (El) 


SfZ[F,0] = 






= - Dct> 


dX0(X)(5F(V)0(A:)J exp (-A[ 0,0]) 

(5 r , . 5 


dX 


< B 


1 

'2 


6J{X) 

dX 6 F{V)Gn{X,Y) 


6F{V) 


' B 


Y=X 


6J{X) 
Z[F,0]. 


Z[F,J] 


j=o 


(A.8) 


Here SF (V) means arbitrary variations of the coefficients of the operator, Sa^^{X), 
6b^{X), 6c{X), and the double arrow implies symmetric action of two first-order deriva¬ 
tives of -F(V) on both arguments of Gn{X, Y) similar to Ea. ()2.2|l 


' B 


Y=X 


dX 5 F{V)Gn{X,Y) 

5a^^{X) - 2 6b^iX) + 5c{X) G 7 v(X, Y) 

= Tr iF(V)G'jv. 


= / dX 

Jb 


Y=X 


(A.9) 
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(A.IO) 


Therefore, from one gets 

Z[F,0] = ~Tt SF(V)Gjv = -^SlnDetjyF = ^<^lnDetGjv 

This expression in the variational form justifies the representation of the prefactor in 
(EHl) in terms of the functional determinant of the Neumann Green’s function (or the 
inverse of the determinant of F(V) subject to Neumann boundary conditions). Si¬ 
multaneously, it gives the variational definition of these functional determinants which 
specifies how the Neumann boundary conditions enter them and how the action of 
differential operators should be understood in the sense of integration by parts. The 
derivation of the gaussian integral ijHB with hxed helds at the boundary can be done 
along the same lines ^HI and leads to a similar variational definition with Gn replacing 
Gjy. 

B. Heat kernel of the square-root type operator 


The heat kernel (j4.dd|l obviously satishes the following Dirichlet problem IC( r j x, ^): 

(—d^ — D + m^)IC(rlx,:i/) = 0, (B.l) 

K(0lx,p) = S(x-p). (B.2) 

Similarly to its solution is given by the ”brane-to-bulk propagator” 

K{t\x,p) =-G\^^\x,t\p,t') W=-dT', (B.3) 

where r | y, r') is the (d + l)-dimensional Dirichlet Green’s function of the 

operator qG+i) = (9^ -|- □ on half-space r > 0 

(m^ - {x,t\p,t) = 6{t-t)6{x,p), (B.4) 

G^;^^\x,0\p,t') = 0. (B.5) 

By the method of images one can construct it in terms of G^'^^\x,t\p,t') - the 
Green’s function in full space without boundary, 

p,t) = G^‘^^\x,t\ p,t)- G^'^^\x,t\ p, -r). (B.6) 

For the massive case it reads 


1 2 

G^^^\x,t\p,t') = Z = m^/\x-p\‘^ + {t (B.7) 

Substituting (|BSl)-(jB21) into (IB.3jl and using recurrent relation between the modihed 
Bessel functions of different orders one obtains (irail . 
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